function x = soltrsup(A,b)

	if det(A) != 0
		N = rows(A);
		i = N;
		while (i >= 1)
			s = 0;
			if (i != N)
				for j = [i+1:N]
					s = s + A(i,j)*x(j);
				endfor           
			endif
			x(i) = (b(i)-s) / A(i,i);
			i = i-1;
		endwhile
	else
		fprintf('ERROR, la matriz ingresada es singular!!, pruebe con otra\n')
    endif
    
endfunction

